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Abstract 

We calculate joint moments of the characteristic polynomial of a random unitary 
matrix from the circular unitary ensemble and its derivative in the case that the power 
in the moments is an odd positive integer. The calculations are carried out for finite 
matrix size and in the limit as the size of the matrices goes to infinity. The latter 
asymptotic calculation allows us to prove a long-standing conjecture from random 
matrix theory. 



1 Introduction 

There is a deep, and still only partially-understood, relationship between analytic number 
theory, and the theory of random matrices. This connection goes back to Montgomery 
[1] who conjectured that statistical properties of non-trivial zeros of the Riemann zeta 
function could be predicted by studying the large asymptotics of correlation functions of 
eigenvalues oi N xN random unitary matrices. This conjecture is supported by theoretical 
[DISlElll], heuristic [3 El [7] and numerical [SIEKIO] evidence. 

The full power of the conjectured relationships between random matrix theory and 
number theory is found in the study of moments of the Riemann zeta function. Using 
random matrix theory, mathematicians have been able to make predictions for moments 
of various kinds, where no conjectures or guesses existed before pT| [T2 l [T3 l [T^ [T5 l [T6 l [T7] . 

A number of review articles have appeared, such as |18l 1191 \2U[ I21j. which summarise 
the main developments that have occurred over the past few years. 

The main object of our study will be the A^-dimensional circular unitary ensemble 
(CUE) of random matrix theory. This is the probability space consisting of the set U(iV) 
of N X N unitary matrices, equipped with normalised Haar measure, fi^^'^^ . 

For a matrix U G U(A^) we denote the characteristic polynomial by 

N 

Zu{9):=l[{l-e'('"-'^), (1.1) 

n=l 

where e'^^, . . . , e^^^ are the eigenvalues of U. 



Define 



Vu{e) := exp(iiV^ - iEli y) Zu{0). (1.2) 



Then VuiO) is real- valued for 9 € [0, 2tt). 

In recent years there has been interest in the joint moments of tlie distribution of Vjj 
and its derivative. Define, for h > —1/2 and k > h — 1/2, 

FN(.h,k):= [ |yt;(0)|2'=-2^^T4(0)|2'^d^«-^ (1.3) 

JV{N) 

and the limiting values 



When h = 0, the moments ()1.3p of V^y are precisely the same as the moments of the 
characteristic polynomial. Keating and Snaith [11] considered F]sf(0,k), and proved that 

FN{0,k) = [[ j^^.^j^y , (1-5) 

and showed that ()1.5p extends to the region D\t{k} > —1/2 of the complex plane. 
Let 2(t) denote Hardy's function: 

Z(0:=e'''WC(^ + it), (1.6) 

where 

m ■■= 3m {log (vr-*/2r (1 + lit))] , (1.7) 

and C{s) and r{s) denote respectively the Riemann zeta function and the Euler gamma 
function. It follows from the functional equation for (^(s) that 2{t) is real for t G M, and 
it is apparent that |C(| + it)| = |^(t)|, so Z{t) is to the Riemann zeta function as Vu{9) 
is to the characteristic polynomial of a random unitary matrix. A series of conjectures 
due to Hall [22j, Conrey and Ghosh [23j, and Hughes [24J has culminated in the following 
prediction for joint moments of 2{t) and its derivative: 

^ r \Z{t)\'^^-^''\Z'{t)\'^^ dt ^ F{h,k)A{k) {logTf+^^' , asT^oo, (1.8) 
T Jo 

where 

p prime ^ ^ e=0 ^ \ J / 

Hughes |24j used random matrix theory to calculate F{h, k) for h = 1,2,3 and Dehaye 
|25t I26j has derived formulae for F{h, k) for all /i € N in terms of sums over partitions (see 
section [2] below for a precise statement). Using their results, the values -F(l,l) = 1/12, 
i^(l,2) = 1/720 and F{2,2) = 1/6720 can be calculated. The corresponding moments 
(jl.Sp for Hardy's function have been calculated by Ingham [2^ and Conrey |28J, and give 
complete agreement for these values of k and h. 

Conrey and Ghosh [29J have also proved (assuming the Riemann hypothesis) that 

1 £ \z{t)z'{t)\ dt ~ ^ (iogr)2 . (1.10) 
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This is proved by relating the joint moment to a discrete second moment of the Riemann 
zeta function at its successive extrema on the critical hue, which had been calculated 
earlier in [30]. The latter result was proved by an integration against the logarithmic 
derivative of a function with zeros at the locations of maxima of + ii)| and which 
could be well-approximated by a Dirichlet series. The numerical constants in (jl.lOp arise 
as values of residues at poles in the relevant contour integral. 

The asymptotic (jl.lOp naturally leads to the conjecture [231 page 110] that 

F(i,i) = 5^. (1.11) 

However, most attention on the problem of calculating moments p.3p has focussed on 
integer values of h. In this work we will take the first steps beyond integer values of h, 
by studying F]\f{h,k) for half-integer values of h. In particular we will supply a proof of 

dm]). 



2 Main results 

In order to put our results into context, we first recall a result of Dehaye |25[ I26j . To do 
this it will be necessary to fix some notations regarding combinatorics of partitions. 

We recall that a partition is a finite sequence A = (Ai, . . . , Aj) with Ai ^ A2 ^ • • • Aj. 
j is the number of parts of A, which is also denoted by ^(A). The sum of the parts of A is 
denoted by |A| = Ai -|- • • • -|- Xj. For n G N we write A h n if |A| = n, and A hfc n if |A| = n 
with £{X) ^ k: a partition of n into not more than k parts. 

The generalised Pochhammer symbol \b]^^'^ is defined for a partition A, a parameter 
£7 > and 6 G M by 

i''ii"'-nn(''+^--i-^)- (") 

i=i j=i \ / 

We will most often be taking the parameter cr = 1, so that we define the special notation 

[&]a:=[6]?^ (2.2) 
In terms of the usual (rising) Pochhammer symbol (•)., we have 

1{X) 

[b]x = \{{b-i + lK. (2.3) 

A partition can be represented graphically by a Ferrers diagram (see figured]), in which 
parts of a partition are represented by a vertical arrangement of boxes aligned at the left- 
hand side. For each box □ in the Ferrers diagram, the arm-length a(n) is the number of 
boxes strictly to the right of □, and the leg-length g(^) is the number of boxes strictly 
below □. The hook-length of □ is a(n) -|- (?(□) -|- 1: the number of boxes to the right 
and below, with the box itself counted exactly once. For the partition A = (4, 3, 1, 1), the 
hook-lengths are indicated in figure [TJ The product of all hook-lengths will be called the 
hook-length of the partition, and denoted hx. For example, for A = (4,3,1,1) we find 
hx = 1680. 
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Figure 1: The Ferrers diagram (left) and the hook-lengths (right) for the partition A 
(4,3,1,1). 



The Ferrers diagram for a partition can be used to define the transpose partition, by 
reflection of the diagram about the main diagonal. So for A = (4,3,1,1), the transpose 
partition is A""" = (4, 2, 2, 1). Clearly the length ^(A""") of a transpose partition is equal to 
the size of the largest part of A, and lA"""! = |A|. It is also straightforward to see that 

Mat = (-1)1^1 [-fe]A and h^T = hx- (2.4) 

We define two quantities in terms of a sum over partitions. Let € N and p € 
No :=NU {0}. Then 



and 



We observe that 



CNip,k) r^C{p,k)NP asAf^oo. (2.7) 

Related quantities appear in the work of Dehaye (see the comment following theorem 12.1 
below) . 

The main result of Dehaye \25\ [26] relevant to our work is the following: 
Theorem 2.1. For h,k eN with k> h-\, 

FNiKk) = i^F^(0,fc)5^^^^^^-^C^(p,fe). (2.8) 



Moreover, for fixed /i € N, Dehaye has shown that the equation (|2.8p extends mero- 
morphically in k to the region lHe{A;} > h — ^ of the complex plane. 

We have re-written Dehaye's result using our notation. In fact he considers quantities 
related to Cn and C (see for example equations (10) and (11) of [26]), which are defined 
similarly to ()2.5p and ()2.6p . but without the restriction on the number of parts of A in the 
summation. The presence of the factor [k]x mean that his and our quantities coincide for 
p ^2k (but could be different for p > 2k). Since the sum in (|2.8p goes only up to 2h, this 
difference is not pertinent in theorem 12. li 

The main result of our work is the following, which gives an explicit formula for -F/v(/i, k) 
for half-integer h: 
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Theorem 2.2. Let h = {2m - l)/2 forme'M and letkeN with k> h-\. Then 

L p=i e=i ^ ^ ^ 

- E (2.) 

p=2/i+l J 

The paper is structured as follows: In section [3]we write down an integral representation 
for FN{h,k) involving an integration over a real parameter and a multi-dimensional 
integral of size A^. In section |3] we evaluate the multi-dimensional integral in closed form, 
and in section O we calculate asymptotics of some integrals related to the ("-integral. We 
give the proof of theorem 12.21 in section [6l and we use the theorem to give evaluations of 
some moments (including (jl.lip l in section [71 



3 An integral representation 

3.1 Notation and properties of Vandermonde determinants 

Let us fix some notation which will be used throughout the remainder of the paper. Let 
X = (xi, . . . ,xn)- We will denote in multiple integrals, 

dx = dxi • • • d^AT. (3.1) 

We shall also denote by 

A(x):= J] (xk-x,), (3.2) 

the Vandermonde determinant. It seems prudent at this stage to note a few properties of 
A(-), that we will make use of later [3lj- 

First of all, note that the Vandermonde determinant is a matrix determinant. We have 
A(x) = det (xt^) (3.3) 

V * Ji,j=l,...,N 

= Y: sign(a)xf)-^..4^)-\ (3.4) 

where Sn is the symmetric group on elements. We see, therefore, that A(x) is a 
homogeneous polynomial in the variables xi, . . . , x^r, of degree N{N — l)/2. 

If {Pji')}j=i,-,N is a set of monic polynomials with the degree oi pj being j — 1, then 
an alternative expression for A(x) is 

A(x) = det (Pi(a;i)),j=i,...,7v ' (3-5) 
which may be proved by applying elementary column operations to the representation 



In order to justify the convergence of certain integrals, we shall employ the following 
crude bound. 
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Lemma 3.1. For x € 



|A(X)| ^ N\ ((1 + x2)V2 . . . (1 + 4)1/2) 1 _ (3_g) 

Proof. For a permutation cr, signo" = ±1. So, by (j3.4p we can bound 
|A(x)l < \xi\''^^^-^---\xnV^'')-^ 

< ( 5^ l) (1 + x?)(^-i)/2 . . . (1 + x^)(^-i)/2, (3.7) 

and use the fact that Sn has order A^!. □ 
3.2 An integral representation for FN{h,k) 

To evaluate the averages Ffyf(h, k) we express this quantity as a multi-dimensional integral: 
Proposition 3.2. Lei n G No, and define 

(_l)n d"- f e \ 

A-„(.,0:=>^^(^ . (3.8) 



T/zen if 2h € Nq anc? k > h - \, 



J — i J 

(3.9) 

In order to prove proposition 13.21 let us collect a few auxiliary results. 
Lemma 3.3. For x € M and e > 0, we have 

/oo 
K„(£,C)e'"'^dC= |xre-^l^l, (3.10) 
-oo 

and the integral converges uniformly in x and e ^ Eq > 0. 

Proof. By calculating residues, or the Fourier inversion theorem, we know that for e > 0, 

ei^<dC = e-"l^l, (3.11) 



uniformly for e ^ Eq > 0. To justify differentiation under the integral, we note that 

^e^ + c\(n+m (3-12) 
for n ^ 1, so that uniform convergence of the resulting integrals is assured. □ 
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To compute averages over U(A^), the most useful tool available is Weyl's integration 
formula. For any function f{U) of a unitary matrix, which depends only on the N eigen- 
values e'^i, . . . ,e*^^, 



/(C/)d^Haar^ 
V{N) (^^J 



l<:j<k<:N 



(3.13) 

By following the substitutions made in [11], we can express F{h,k) as an integral over 



Proposition 3.4. It follows from ()3.13p that 

nN^+2kN-2h foo foo N 



FN{h,k) 



Proof. Differentiating ()1.2p . we get 



: |xi + ••• + X7vp''A(x)2dx. (3.14) 



So, 



Furthermore, 



VIj{0) = Vu{0) 



z'uie) _ d 



2 ^ Zu{0)J' 



ZuiO) d.l°s(^-(^)) 
^J!, iei(gn-e) 
~ 2^ 1 _ ei(fn-e) ' 

n=l 



from (jl.ip . We then get 



N 



using the fact that 



n=l 

^-^%(o)|;cot(^). 



e-ie„ _ 1 2 ^° V 2 / 2 ' 



Substituting ()3.18p into ()1.3p . we get 

nh,k) = ^[ \Vu{0)\ 



\2k 



U(7V) 



2h 



Haar 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



Following |llj . we write the integral over U(A^) as a multiple-integral using Weyl's identity 
(|3.13p . and make the substitutions Xj = cot {6 j/ 2) therein, leading to p.l4p . □ 
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Proof of proposition [372l By lemma [331 for any e > 0, 

/CO /"OO 1 

• 7_ n (TT^iyra + ■ ■ ■ + -^1'' e-l-+-+-l A(x)2 dx. (3.21) 

By lemma 13.11 and equation ()3.12p we see that the left-hand side of ()3.2ip is absolutely 
integrable, so by the Fubini-Tonelli theorem, we can move the (^-integral to the outside. 

Considering the right-hand side of ()3.2ip . we estimate 

= N^^\l + xlf---{l + xj^f. (3.22) 
Combining this with lemma 13.11 we get 

TV ^1 

so that \i k > h — 1/2, the right-hand side of ()3.2ip is uniformly convergent in e, and we 
may pass the limit e J, under the integration. The resulting equality is (|3.9p . □ 



4 Evaluation of the mult i- dimensional integrals 

The main calculation in this paper is an evaluation of the integral 



/OO /"OO 1^"^ * 

.^••7_^n (i + ^2jiv+. ^W'dx. (4.1) 

To this end, we will first derive an equivalent representation for H{k, C,) which we will 
be able to evaluate, for integer fe, in terms of Laguerre polynomials, and multivariable 
hypergeometric functions. 

Proposition 4.1. Let ( > be fixed, and A; € C with 9^e{A;} > — ^. Then 

N N-1 /-oo /-oo ^ 

(4.2) 

Proof. Let 

roo poo ^ 

L{k,C):= •••/ l[{y,+2C)''yyy^A{yfdy. (4.3) 
Jo Jo jj-^ 
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By ()3.5p we may write, for fixed C, 

= ^ sign((7)(yi + 2Cr^'^'' ■■■{yN + 2Cr^^^-\ (4.4) 

Using the expression ()3.3p for the second of the two Vandermonde factors, we can express 
L{k, C) as a sum of products of integrals. We get 

N 

L{kX)= X] sign(cr)sign(r) (4.5) 

o-,tG5jv j = l 

where 

poo 

I^,.:= (2/ + 2C)^+V+'e-^dy. (4.6) 
Jo 

Let us now consider the integral 
We exploit the homogeneity of A(x) to write 

{Xm ) 



,^ , , n — i-T 



j=l ^ 0' Is^jKmfiN 



TT -^i TT 

li (1 + + ix„,) „ (1 - ixj)a - ix„) 



n n (— —) 

\l + ixm l + ixnj V 1 — ia^m 1 — ix, / 



E 



sign((T) sign(T) 



^ (1 + ixi)-(i)-i(l - ixi)-(i)-i • • • (1 + ix^)-W-i(l - ixw)-W-i ■ 

(4.8) 

Substituting this into ()4.7p we get 



• n (1 + ^2^N+k ^ W = sign(fT)sign(T) J]Xo)-i,rOVi, (4.9) 



where 



j = l 1' cr,rGSjv J = l 



^'"'^ '~ i-oo (1 + ix)i+'=+-(l - \xY+^+^^ ^^-^^^ 

The integral (|4.1Up can be evaluated in terms of a confluent hypergeometric0 function — 
from |321 formula 3.384.9] we find that for a, /? G N and for C > Oj 

e'*^^ TT ria + B - I') A 

/ 7 — dx= ^ t,, e~^iFi(l-a;2-a-/3;2C) 

J_^{l + ix)^{l-\xY 2"+/'-2 r(a)r(/3) i iv . s) 



^For the definition of i_Fi, we refer the reader to section [42 
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— however, our main concern is to show that J^^u is equal to up to a constant, which 
we do next. 

We, temporarily, assume that d\z{k] > 0. Since 

1 1 /"°° 

7 ^TT^ = ^, ; ^ e-^^-^^^V+'^dn, (4.11) 

we may write 

1 poo poo piC^— C 



We bound 



and 



1 



(1 + ix) 



l+k+u 



g(7r/2) 3m{fc} 

^ (1 +2;2)(i+'^+SKKfc})/2' ^^-^^^ 



^-(l-\x)u^k+^l 



so that the double integral (|4.12p is absolutely convergent and we may reverse the order 
of integration. Therefore, we have 



r{l + k + fi) Jo J_^ (1 + ix)i+*^+- 



2tt 



r{i + k + fi)r{i + k + u) 



22fc+M+'^r(i + k + /i)r(i + k + u) 

TT 

22k+t,+upi^l ^f.^ + k + u) 



poo 

Jo 

poo 

/ (2/ + 2C)''+V+'^e-J'dy 
Jo 



(4.15) 



using Laplace's formula (equation 3.382.6 of 



s-l -p 

dx = ^ ^ , . , p > 0, (4.16) 



27r (1 + ix)^ r{s) 

to pass from the first to the second line of ()4.15p . 

The relationship (1^5]) . together with (gSD and (g^]) proves (02]) for lHe{A;} > 0. To 
complete the proof we note that both sides of ()4.2p may be continued as analytic functions 
oiktodlc{k}> □ 

For integer values of k we are able to give two direct evaluations of the integral in the 
right hand side of ()4.2p . The first one uses Laguerre polynomials, and the second uses a 
hypergeometric function of matrix argument. 



4.1 Evaluation in terms of Laguerre polynomials 

We 

by 



We recall that the classical Laguerre polynomials are defined for a parameter a > —1 



LH(t):=^il(t«+ne-*). (4.17) 
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An explicit formula for L^f^ is given by 

r(«)rti-V" r(n + Q + i) {-ty 

- j^, rO- + a + l)(n-,)!^- ^^-^'^ 
In terms of hypergeometric functions there is the following expression |331 formula 22.5.54]: 

^'^^W = + (4.19) 

Applying Kummer's transformation [33l formula 13.1.27] to (j4.19p leads to an alternative 
expression, 

Lio.) _^t^ r{a + j + n) i-ty 

We denote by W{gi, . . . , gn){x) the Wronskian of the n functions gi, . . . , gn, evaluated 

( 91 
9'i 



>V(ffi,---,5n) := det 

We have, 

Proposition 4.2. For A; € N and ( G 



I (n 



92 


9n \ 


9'2 ■ 


9'n 


(n-1) 
92 





(4.21) 



H[K C) = (-l)'=('=-^)/^|?0e-^ICI>V(LSJ), lSJ) „ . . . , 4Vi)(-2|CI). (4.22) 

Proof. The integral on the right-hand side of ()4.2p is the averaged moment of the charac- 
teristic polynomial of random matrices from the Laguerre unitary ensemble. Such averages 
were considered by Brezin and Hikami, who showed [341 page 114] that 

roo roc ^ N^Y\^S^C- 

/ •••/ ll(*-yj)'y^e-y^A{y)'dy= ^^if ' W{pn , ■ ■ ■ ,PN-,k-im, (4.23) 
Jo Jo Ue=o^- 

where pj = {—ly jlL^f'\ and 

roo 



= jl^k\(^^l^ ) =jl{j + k)l (4.24) 

(the last line being a classical result of Laguerre polynomials [35]). We have 

fe— 1 

W{pN, . . . , PN+k^i) = (-l)^^+^(^-i)/2 J] (AT + ,)! w(lS) , . . . , L^^l,_,) 

j=0 
N-1 

^ (_l)fcAr+.(.-i)/2 -Q (^k + jy.WiL^P,...,L^;;l,_,) (4.25) 

j=Af-A; 
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and 



niV-l „ N-1 N-1 
_J=0 

Y[£=0 ^- j=k j=0 



k~l 

V. 

N-k-1 N-1 



= n n(-?'+^)' (4.26) 

j=0 j=0 

Putting (10511 and ()i:26]) into (I03]l . and setting t = -2C gives 
+ 2C)'yK^^A(y)2dy 

i=i 

N-1 

= ^.l)H>^-mNiY[U + k)l'W{L^i\L^il,,...,L^^^^^^^ (4.27) 

j=0 

Together with proposition 14.11 this proves ()4.22p for ^ > 0. For general (" S M, we note 
that the function H{k, () defined by (|4.ip is an even continuous function of C, so that we 
may replace by \C\ wherever it occurs. □ 

For k = 1, the evaluation (|4.22p reduces to 

H{1,0 = ^^e-^l^lL«(-2|Cl), (4.28) 

so that the integral ()4.ip is proportional to a single Laguerre polynomial. Although we 
were not able to find a use for (j4.22p for general k, we present in appendix[X]an elementary 
proof of (jl.lip . based on ()4.28p . and properties of Laguerre polynomials. 

We are also able to write the Wronskian appearing in ()4.22p as a Hankel determinant 
without derivatives, which may be of independent interest. 

(k) 

Proposition 4.3. // L]y denotes a Laguerre polynomial, then 

W{lP,. . . = det ■ (4-29) 

Proof. We make repeated use of the identity 

^^Lt\t) = -Lt\'\t), (4.30) 



to get 



w(LS;),...,L(;^,_,)(t) = det(^L!;i,(t)) 



i,j=0,...,k-l 

det((-l)^Li;:£/t)^ 



i,j=0,...,k-l 



= (-l)-det(4:e,w),^^.^„^ (4.31) 
where w = A;/2 if A: is even and (A; — l)/2 if /c is odd. By application of the identity 

Ll^-'\t) + Li%{t) = Lt\t) (4.32) 
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via row operations on the matrix in ()4.3ip . we get 

W(i!f^•••^4Vl)W = (-irdet (Lf-'J^{t)), (4.33) 

Finally, switching rows according to j i— )• k — l—j gives uj transpositions, leading to ()4.29p . 

□ 

4.2 Evaluation in terms of hypergeometric functions of matrix argument 

The (single-variable) hypergeometric function pFq is defined formally [36j by the series 



.F,(ai, . . . , ap; 5i, . . . 6,; x) := f; Ml^TT^ ^^.34) 

^ (^i)i • • • (feg)j 



where (■). is the rising Pochhammer symbol. The parameters oi, . . . , Op, bi, . . . ,bq can be 
arbitrary complex numbers, however if 6j E Z\N then the series (|4.34p becomes undefined 
unless there is a coresponding parameter Oj' G Z with bi < Oj/ ^ 0, in which case we adopt 
the convention that the series (j4.34p terminates after Oj/ terms. 

The hypergeometric functions of a matrix argument provide a multi-variable generali- 
sation of (|4.34p . They have been studied in [371 ESI |39] , and have been found to occur in the 
context of random matrix theory in the statistics of extreme eigenvalues [40 ^ HTj H2| HH] 
and moments of characteristic polynomials off the critical line [44 [ I45| . amongst other 
places. To generalise ()4.34p . the sum over integers in is replaced by a sum over partitions, 
the Pochhammer symbols are replaced by the generalised Pochhammer symbols defined 
in section [21 and the univariate monomials are replaced by Jack polynomials (see [l6] 
or [m chapter 12]). 

Let o" > be a parameter, and X be an x matrix with eigenvalues xi, . . . ,xn. 
Then 

pP'fXai, ...ap;bi,...,bg;X) := pF^f^ai, . . . Op; 6i, . . . , 5^; xi, . . . , xn) (4.35) 

where the Jack polynomials are [46j homogeneous symmetric polynomial eigenfunc- 
tions of the partial differential operator 

2 ^ ' dx^ ^ Xi- Xj dxi 

i=l ' iytj ■> 



normalised so that 



cr(i,...,i 



a 



(4.38) 

A 



(this is different to the normalisation used in |46j ) where c\{a) and c\{(t) generalise the 
hook lengths defined in section [21 their definitions are 

cxi^) ■■= n (^(°) + 1 + ^^(°)) ' (4-39) 

□ SA 
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and 

c',{a):=Yl{i{D) + a{l + 9m). (4.40) 

□ga 

For the case a = 1 the polynomials are proportional to Schur polynomials. In that 
case, the denominator in (|4.38|) is c\{l)c'^{l) = h\. 

The convergence of the series in (j4.34p and (|4.36|) depends in general on the parameters 
ai, . . . , Op, 61, . . . , 6q. If p — g ^ 1 then the radii of convergence in (|4.34p and ()4.36p are 
at least 1. If G Z \ N for any i then the series terminate, and are thus defined for all 
values of the arguments (in which case the resulting functions are polynomials) . 

Hypergeometric functions of matrix argument enjoy a reflection property |48[ page 812] 
that we will make use of. Let a € N and b > a. Then, 

^F["\-a; -b;xi, . . . ,x„) = 



" r(6-a + i^ + l 



n-1 -1 -1 



f=i r(6 + ^ + l) V a XI X, I 



This may be viewed as a generalisation of the single-variable identity: 

r{b-a + l) , / , , , 



iFi (-a; -b; x) = ^(^^1) (^-«' ^ + ^ " V j " ^^'^^^ 

In both (j4.4ip and (j4.42p we emphasise that since the parameter a is a negative integer, 
the hypergeometric series are actually finite, and the left-hand side series terminate before 
the denominators in the summands in ()4.34p and ()4.36p become zero. 

The main result of this section is the following: 
Proposition 4.4. For k G'N and ( ^M., 

^(^' = 2(/+l^i)iV -^^(0' ^)«"'^'^' i^S'^(-^; -2^; 2|C|, . . . , 2|C|). (4.43) 

The proof of proposition 14.41 is based on the following evaluation of the integral L{k, (^): 
Proposition 4.5. Let L denote the integral in (|4.3p . and let k £ N and C > 0. Then 

L{k,C)={ nr(j + 2)r(fc + l+j)j(2C)'=Vi'^(-^;iV + fc;;^,...,^). (4.44) 



Proof. Forrester and Keating \4b\ equation (3.2)] have provecH the following integral, 
valid for 9^e{a} > -1, d\e{b} > -1 and either \t\ > 1 or 2^ G N: 



1 fi fi ^ 



j2^A'^p(l/7) ^_2^, ^,(jv - 1) + o + 1; 2',{N -l) + o + 6 + 2;i,,,,i), (445) 



^There is a typo in equation (3.2) of [JS] in the second parameter of the hypergeometric function. We 
have given here a corrected formula. 
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where 

We let 2^ = A; G N and a = k, b = L, j = 1, t = — 2(^/L, and make the changes of 
variables xj = yj/L. The integral on the left-hand side of ()4.45p becomes 



(-1) 



kN rL f-L N 



Since (1 — yj/L) ^ —> e~^J as L — > oo, and 



/ ny,'(l-|)^y.+2C)'A(y)'dy. (4.47) 



l-f)^l[o,L](y)^e-^^ (4.48) 



for Uj ^ 0, we get, by the dominated convergence theorem, 



L(A;,C)= lim (2C) 



kN 



X L^'+^'^'SNik + 1, L + 1, 1) ^F^P (-k,N + k;2N + k + L;^,...,^) , (4.49) 



where the hypergeometric function of a matrix argument is a multivariate polynomial since 
/c S N. By Stirling's formula and ()4.46p it follows that 



hm L'^'+'^^SNik + 1, L + 1, 1) = TT r{j + 2)r{k + j + l). (4.50) 

L— !>oo -■■ 

j=0 

For a partition A, it is clear that [2A^ + k + L]x L''^' as L — t- oo, and the fact that C^^ 
is homogeneous of degree |A| suffices to conclude from (j4.36p that 

lin.^,F?l^-k,N + k;2N + k + L;^,...,^)=,F^^^I^-k,N + k;;^^ 

(4.51) 

This completes the proof. □ 



Similar integrals have been evaluated in [48j by a different method. 

Proof of proposition \4.4[ Combining the result of the proposition 14.51 with the reflection 
formula (|4.4ip and equation (|4.2p leads to (|4.43p . for C > 0. To pass to the case ( &M.we 
use the same argument as at the end of the proof of proposition 14.21 □ 



5 Some integrals involving Kn{s,C) 

In this section we will consider for integer values of n and p, integrals of the form 

/ if„(£,C)e-^fC^dC (5.1) 
Jo 

in the asymptotic regime e J, 0. 

The function Kn is defined in (|3.8p as a partial derivative with respect to e of a function 
of e and It will be convenient to obtain alternative formulas for Kn with derivatives 
with respect to (. 
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Lemma 5.1. Let Kn{e,() be defined by (|3.8p . If n = 2m is even then 
whereas if n = 2m — 1 is odd then 



vr (9C W + C 

Proof. We begin by observing that 

(-1)" (9" / e 



Kn{e,C) 



(-1)" 5" / 1 



+ 



1 



27r ^e*^ Ve + iC ^ - 
n! / 1 1 



27r Vi"+nC-i£; 



n+l 



-i)"+i(C + ie)' 



n+l 



If n = 2m — 1 is odd, then 



.n+l 



27ri2™ V(C-i£ 



1 1 

+ 



(c + i^: 



+ 



1 



27r 5C™ VC-ie C + ie 
(-1)™+! 5" / C 



On the other hand, if n = 2m is even, then from (|5.4p we get 

a" / 1 1 



27ri 5C VC-is C + is 



Lemma 5.2. Lei L^"^ denote the Laguerre polynomials, defined by ()4.17p . //p 

a" 



9C' 



For p ^ n, we have 



Proof. The case p ^ n is slightly the simpler and we consider it first. We have 



an , 

n!C^-"e-^«4P-")(iVC)- 



16 



If p ^ n, we still have 



(5.10) 



but now we write 



dc 



gn-p 



g^n-p 

At this point we use ()4.20p to write, for a € No, 



i=NC 



(5.11) 



^f.)(0 ^,y^ (a+P + ^-l)!(-0^ 
^'^1 2^ (a + ^)!p! ' 



£=0 



(5.12) 



where the power series converges for all values of ^. Thereby we get an expression for 
e~^L^\^) which we may legitimately differentiate term-by-term to get 



gn-p 

g^n-p 



E 

l=n—p 



{p + l-l)\{-lY £! 

i\ {i-n + p)\ 

(m + n-1)! (-^)™ 



m=0 



(m + n — p)!p! m! ' 



via i = m + n — p, 



= (_l)n-Pe-«4"-?')(e), (5.13) 
using ()5.12p once more. Substitution of ()5.13p into (jS.lip completes the proof. □ 



For the case p > n in ()5.ip . the asymptotic evaluations are given by the following 
proposition. 



Proposition 5.3. Let p ^ n + 1. Then, as e ^0, we have 



f "'^^:""'^ -l)("+^/^+o(l), nodd, 



o(l), 



(5.14) 



n even. 



Proof. We use the representations of Kn{£,C) derived in lemma ETH and then integrate 
by parts, and insert the derivative formute from lemma [5^21 

If n = 2m — 1 is odd then 



i^„(£,C)e-^kPdC 



IT 



5C" w + c 



vr 



£2 + g(^n 



e-^^CdC 



^ 1)™ /"°° ^p—n+l 



-n\ 



TT 



+ 



-^^L(f-")(iVC)dC. (5.15) 
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Since p—n+1 ^ 2 and since | / (e^ + ) | ^ 1 the integral in (|5.15p is uniformly convergent 
in e, so we can pass the limit under the integral and get 

lim / Knie,C)e-^'^edC = - —nl C^-^-V^^L^f"") (iVC) dC- (5.16) 

£^0 Jo ^ Jo 

At this point, we insert the expansion (|4.18|) for the Laguerre polynomial, to find 

\m ^ / ^ \ f \^ /"oo 



£=0 

(5.17) 

(_l)»^n!(p-n-l)! 

- ^ j;^n ' (5-18) 

using that the sum in ()5.17p evaluates to 1, a fact which is proved in appendix IB. 1[ 
In the case that n = 2m is even we again use lemma 15.11 and lemma 15.21 to find that 



f 1 \m+n roo 

= n! / ^^C^-e-^C4P-n)(^^)dC 

71" Jo ^ + C 

= o(l), (5.19) 

as e J, 0, using lemma 15.41 below. □ 

In the proof of propositon 15.31 we used the following standard result (quoted without 
proof): 

Lemma 5.4. If f is a bounded function with f{x) — > /o as x J, 0, then 

/°°^2^/(^)dx = ^ + o(l), (5.20) 

as e 10. 

Proposition 5.5. If p ^ n and n = 2m is even, then 

/ A:„(£,C)e-^^C^dC= ^ ^ , , (-A^)"-^ + o(l), (5.21) 
Jo 2 {n-py. 

as e 10. 

Proof. Applying lemma 15.11 and lemma 15.21 we get 

' / 1 \m,+n /'oo on 

E:„(e,C)e-^k^dC = / ^^i^(e-^^OdC 



vr Jo + C 9C" 



p!(_iV)n-P / ^^e-^^4"-P)(ArC)dC. (5.22) 
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Applying lemma [531 we find in the limit e J, 0, 

/ K,(e,C)e-^«C^dC= p!(-7V)"-^'-4-?')(0) + o(l) 

Jo vr 2 

(_-\\m+n I 

using the explicit representation (|4.18p for the Laguerre polynomial. □ 

For odd n ^ p, we have the slightly more subtle result: 

Proposition 5.6. Let /(C) :=X]p'=o fpC^' with P ^ n, where n is now odd. Then provided 
that 

P 



V(M^ = 0, (5.24) 

then 



// condition (j5.24p does nof hold, then the limit in (|5.25|) diverges. 



Proof. Following the arguments above, we find for n = 2m — 1 and p ^ n 
Since 



i^„(e,C)e-^^C^dC = ^-^p!(-Ar)"-P / _^e-^';L("-^')(iVC) dC. (5.26) 



we have that 



4"-p)(ivc)dc= ! / ^^dc (5.^ 



e2 + c2 p ^ ^ p y 7g e2 + c2 







passing the limit e J, in a similar way as in the proof of the first part of proposition 15.31 
Evaluating the integral on the right hand side of (|5.28|) we get 



+ ^^p!(-iV)-pf^( -^]iJ)!+o(l). (5.29) 

It is clear that when condition ()5.24p is satisfied, the contributions coming from the first 
term of the right-hand side of (j5.29p cancel, and we arrive to (j5.25p . 

If condition (j5.24p does not hold, then we need to prove that the integral on the right- 
hand side of ()5.29p diverges as e J, 0. To see this, we integrate by parts to get 

[■CO A AT r-oo 

/ ^T72«"''^dC = -loge + ^ / log(e2 + CV^^dC (5.30) 
and observe that the integral on the right-hand side of (|5.30p is 0(1) as e J, 0. □ 
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6 Moments of characteristic polynomials and their deriva- 
tives 



We have now collected ingredients required to prove theorem 12.21 After giving the proof 
below, we then will consider the asymptotics limit N oo oi large matrix size (proposition 
U below). 



6.1 The finite N case 

Proof of theorem \2.2\ We use proposition 13.21 and proposition 14.41 to write, for /c G N, 

F^(/i,A:) = lim-2^F^(0,fc) / /^2h(£, Oe""^!^! iFS^^(-A:; -2A:; 2|C|, . . . , 2|C|) dC (6.1) 

We shall derive an expansion for the hyper geometric function in terms of the coefficients 
Cn{p-, k) defined by (|2.5p . We begin by using the definition (j4.36p to write 

,F?i-k; -2k- 2ICI, . . . , 2ICI) = 5: [Ltiu'^^''^''|A|V''''^'^ - ^''^^ 

A 

The series in (|6.2p is finite. Indeed, the sum runs over only those partitions A with largest 
part not greater than k (so that the factor [—2k]x in the denominator is never zero) and 
the Jack polynomial c|^'' vanishes if A is a partition with more than parts. This places 
an upper bound of kN on the sum of the parts of A. We use the homogeneity of Jack 
polynomials, and the normalisation (|4.38p to get 

-2k; 2|C|, . . . , 2ICI) = E ^mr^^^^ 

We can index the sum in (|6.3p by transposes of partitions, rather than the partitions 
themselves. This gives 

,F«(-fe; -2k; 2ICI, . . . , 2|C|) = E r^^(2|Cl)'^l 

using (|2.4p . We will group the terms of ()6.4p so that partitions of the same integer p are 
summed together. Since we transposed the partitions, we know that A can have at most 
k parts. This manipulation brings us finally to 

yp{-k; -2k; 2ICI, . . . , 2ICI) = E f E wF^l (-'i^')' 

p=0 VAhfcp ^ A / 

kN 

= ECiv(p,fe)|Cr- (6.5) 

p=0 
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Substituting ()6.5p into ()6.ip we arrive at 




We split the sum into two contributions according to ^ p ^ 2h and 2h < p ^ kN, and 
apply proposition 15.31 to the second sum, and apply proposition 15.61 to the first sum. By 
proposition 13.21 we know a priori that the limit e J, exists in ()6.6p . so condition ()5.24p 
must hold with fp = Cj\f{p, k). This means that in addition to proving (j2.9p . we have also 
proved the combinatorial identity 

2h 

E( p'jc^fe^^^O. (6.7) 

valid for 2h an odd integer with < h ^ k. □ 

We remark that starting from equation ()6.6p (which does not depend on the parity of 
2h), and using propositions 15.31 and \5.b\ we can re-prove Dehaye's result, theorem 12. H for 
2h even and k integer, using our methods. 



6.2 The N ^ oo limit 



In order to pass to the limit A/^ ^ oo in ()2.9p we require an estimate on the size of the 
coefficients C]\[{p,k). This is provided by the following lemma. 

Lemma 6.1. Let N ^ 1 and p ^ 2. Then 



where the implied constant may depend on k, hut is independent of N and p. 



Proof. We use the fact that if A = (Ai, . . . , A^) is a partition of p into not more than k 
parts, then Ai ^ [|J , where [-J is the integer-part function. Then we have, 

fc Ai Ai 
[2fe]A = llll('^k + j-i)^l[{2k+j- 1) 

i=i 

[p/k] 

> n 

_ r{2k+[p/k\) 



i=ij=i 



I) 



r(2k) 



With this inequality, and the trivial estimate 

\[-NU<^{N + kY, 

we can bound C]\f{p, k) as follows: 



\CM(p,k)\=2-' 



Al-fcP 



[k]x[-N]x 



[2kUl 



^ 2P{N + krr{2k) ^ [% 
^ r{2k+ip/k\) ^ hi 

2P{N + kfr{2k) ^ [k]x 



r{2k+\p/k\) ^ hi 



(6.9) 
(6.10) 



(6.11) 
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where the last equahty holds, since [fc]A = if A is a partition with more than k parts. It 
follows from the hook-content formula (see Theorem 7.21.2 of [IS] combined with Propo- 
sition 2.2 of [50]) that 



Ahp 

SO we have proved 



,^ . ,M NP f kY (2k)Pr(2k) 

which furnishes the required estimate. □ 
Lemma |6. II is probably far from optimal, but is sufficient to prove the result following. 

Proposition 6.2. Let h = {2m - l)/2 /or m G N and let k e'N with k> h-\. Then 

2h p 



I p=i i=i ^ ^ / 



where 



Proof. We recall that 



i2h)\ip-2h-l)\Cip,k)\, (6.14) 

p=2h-\-l J 

k 

fj^ r{k + j) 

F{h,k):=Y.rr.^j^^F^{h,k). (6.16) 

We can apply the limit term- by-term in the first sum of ()2.9p . using ()2.7p and the fact 
that 

FN{Q,k) ^ F{id,k)N^\ (6.17) 
as N ^ oo, which was proved in |llj . 

For the second sum, lemma 16.11 gives us 

^""^^''^■ CNiP, k) = O, (^) , (6.18) 

so that the second summand of ()2.9p is bounded independently of by a summable 
function of p. Taken together with ()2.7p . this allows us to apply Tannery's theorem [511 
§49] to prove 

( — Oh — ^]\ °° 
^ NP ^ ^^(P'^)= E iP-'^h-l)\C{p,k), (6.19) 

p=2h+l p=2h+l 

and hence we get (|6.14p . □ 
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p — X boxes 






1 




X — 1 


X 


x + 1 




p — X — 1 


-1 







x-2 





X boxes 



Figure 2: Partitions of p into 2 parts. The contents of each box are the values of j — i 
appearing in the definition of the Pochhammer symbol. 



p — X + 1 


p — X 




p-2x + 2 


p — 2x 




2 


1 


X 


X — 1 




1 





Figure 3: The hooks for the partition X = {p — x,x). The hook length hx is the product 
of the entries in the boxes. 



7 Partition sums and a proof of ( II. Ill) 



In order to give explicit formute for the moment F{h, k) we require closed forms for C{p, k). 
We have been able to find these forms for k = 1 and k = 2: 



Proposition 7.1. Let p € N and C{p,k) be defined by (|2.6p . Then, 

2P 



C{p,l) 



and 



Proof. We recall that 



C{p,2) 



12(2p + 4)!2P 



pi{p + 2)\{p + 3)\(p + 4y: 

[k]x 



(7.1) 
(7.2) 



C{p,k) :=2P J] 



where the summation goes over partitions of p into not more than k parts. In the case 
k = 1, then only the single partiton A = {p) is admitted. In this case, it is easy to see that 
[1]a = p!, [2]a = (p + 1)! and hx = p\. This immediately leads to ()7.ip . 

For k = 2 we require all partitions of p into not more than 2 parts. These partitions 
are of the form A = (p — x,x), where ^ x < (p + l)/2 (see figure [2|). For a partition of 
this form, we have 

[2]a = {p-x + l)lx[ 

and 

{x + 2)l{p - x + 3)\ 



[4]; 
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(7.3) 
(7.4) 



and the hook length hx is given by (see figure [3]) 

x\{p — x + 1)! 



hx 



so that 



Cip,2) = 2P ^ 



0<a;<(p+l)/2 



p — 2x + 1 

12{p -2x + 1)^ 

xl{x + 2)!(p - X + 3)!(p - X + 1) 



(7.5) 
(7.6) 
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We observe that the summand in ()7.6p is invariant under the reflection x p + 1 — x, so 
that the half-range sum can be replaced by one half times the sum from to p+1, giving 

C{p, 2) = 2^6^ ^ 3. + 3)1!^ - X + 1)! ■ ^^-^^ 



x=0 



In appendix IB . 21 the sum in ()7.7p is evaluated, whereupon ()7.2p follows. □ 

Based on the results of proposition I7.lt it is tempting to conjecture that C{p, k) will 
be 2P times a ratio of products of factorials for all /c G N. However a computer-based 
investigation of C{p, 3) has shown that this structure appears to break down when A; = 3. 

When h = 1/2 and k = 1, substituting ()7.ip into equation (|6.14p gives 

2P(p- 2)! 



,2 



5 



47r ' 

as conjectured. 

Using (|7.2p for the case A; = 2 in (j6.14p . we can evaluat^ 

^(i 2) = (15 - 73F3 (1, 1, 1; 3, 6, 7; 8) ) 



(7.8) 



0.00815... (7.9) 



and 



= iMo^ (^^^^^ ^' ^' ^' ^' ^' " 

« 0.000 354... (7.10) 
These last two values do not appear to have been derived or conjectured before. 
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A An elementary evaluation of 1) 

In this appendix we give a completely elementary derivation of (II. lip based on classical 
properties of Laguerre polynomials (see, for example, |35j). 

^We have used the algebraic manipulation package Maple to derive closed forms for the sums in terms 
of hypergeometric functions. 
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We first recall the alternative expression (|5.3p for Ki{e,C)- 



Kiie,C) = z^J-T^]. (A.l) 



Equipped with this, we proceed from the integral representation of proposition 13.21 and 
equation (|4.28p . 

= lim- / Ki(£,C)e-^l«lL!i^(-2|C|)dC. 
Substituting (jA.ip and integrating- by-parts, we get 

= if K« (-20 - <'-i(-2C)) dC, 

where we have used the fact that 

ALS,")(t) = -LK)(t). 
By a standard recurrence for Laguerre polynomials, 

TVLff (-20 - 2LgU(-2C) = 2CLg)_i(-2C), 



so that we have 



Since 



F.(il) = -nmy^ ^e-C^SJ)_,(-2C)dC. 



+ C"' 

we can apply the dominated convergence theorem to pass the limit under the integral, to 
get 



o f'OO 

F^(i,l) = -/ e-^^LgLi(-2C)dC. (A.2) 



At this point, we insert ()4.18p . the explicit representation of L^^\ getting that 

/ ,x . „ \ r)n poo 



N + 2\ 2" 



So, 



n=0 

^ V n + 3 y AT^i+i' 
n=0 



n=0 ^ ^ 
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and 

n=0 ^ ^ 

/ iV + 2 \ (A^ + 2)''+3 
To perform the last sum, we use the estimate o ^ — 7 — Set 



n + 3 J (n + 3)! 



N + 2\ 2" ^ 2" /^Af + 2^ 



n + 3 J Arn+3 - (^ + 3)1 \^ TV 
2" 2 
n! 

which is summable. This ahows us to use Tannery's theorem |5H §49]. Since 

1 ^ + 2 \ 1 
N^oom+^ \ n + 3 J (n + 3)!' 

by Stirhng's approximation, we get 

2 v— r 2*^ 

n=0 ^ ^ 



^2 



5 
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B Some sums 

In this appendix we give evaluations of some finite sums that appear in the main text 

B.l A sum from section O 

In this appendix we evaluate the sum 

P \ I p — n + £ — 1 



(A.5) 



n-i 

for p ^ n + 1. 

By elementary manipulation of binomial coefficients, 

P \ f p — n + £+l\ pl{p — n + £-iy. 



n-£j\ £ J (^n - £)l{p - n + £y.£\{p - n - 1)\ 

p ( P — ^ \ ( n 



p - n + t \ n 

We also use the fact that 



(B.2) 



1 1 

p — n + £ 
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together with ()B.2p to get 



n 







We have used 



/ -x)''dx= r— ^> fora,6€N: (B.5) 

Jo (a + 6 + 1)!' 



a version of Euler's integral. 

B.2 A sum from section [7] 

In this appendix we evaluate the sum 

^x!(x + 2)!(p-x + 3)!(p-a; + l)!' ^ ' 

Let us observe that the numerator in ()B.6p can be written as 

{p~2x + lf = {p-x + l){p-x + 3) + x{x + 2)-{p-x + l){x + 2)~x{p~x + 3). (B.7) 
Therefore, 

^x!(x + 2)!(j9-x + 3)!(j9-x + l)! ^ ' ^ 

P+i ^ 



-j^ . - 



x=0 



x\{x + 2)!(p - x + 2)!(p - x)! ^ (2; - l)!(x + l)!(p - x + 3)!(p - x + 1)! 



' ^ + l)!(p - X + 3)!(p - x)! ^ (a; - + 2)!(p - x + 2)!(p - x + 1)! 



(p + 2)P \^\x + 2l\ X I ^ \ X -1 I \ x + 1 

^■f^ ' \x=0 ^ ^ ^ ^ x=l ^ ^ ^ 

1 (^( p + i\(p + ^\,Y^(p + i\(p + ^ 

(p+l)l(p + 3)l[^^\x + l J\ X )^Z^^{x-l J\x + 2 
2 f p + 2\ f p + 2 



E 



(p + 2)!2 ^ \ X J \ p-x 

1 (^(p + 3 \ ( p + l \ ^ ( P + l\ ( P + 3 

{p+l)l{p + 3)l[^^\ X )\p-x) ^ )\p-x 
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We then use the Vandermonde identity, 



x=0 

to get 
p+i 



{p-2x + iy _ 2 f 2p + 4 

P 

2 / 2p + 4 



^ xl{x + 2)\{p -x + 3)\{p -x + l)\ {p + 2)!2 V P 



{p+l)\{p + 3)\ \ p 

1 1 W 2j5 + 4 



{p + 2)!(p + 3)! V P 



(p + l)!(p + 2)! Vp + 2 P + SJ\ p 

2p + A 



(B.IO) 
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